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Abstract 

In a recent work on the quantization of a massless scalar field in a particular 
colliding plane wave space-time, we computed the vacuum expectation value of the 
stress-energy tensor on the physical state which corresponds to the Minkowski vacuum 
before the collision of the waves. We did such a calculation in a region close to both the 
Killing-Cauchy horizon and the folding singularities that such a space-time contains. 
In the present paper, we give a suitable approximation procedure to compute this 
expectation value, in the conformal coupling case, throughout the causal past of the 
center of the collision. This will allow us to approximately study the evolution of such 
an expectation value from the beginning of the collision until the formation of the 
Killing-Cauchy horizon. We start with a null expectation value before the arrival of 
the waves, which then acquires nonzero values at the beginning of the collision and 
grows unbounded towards the Killing-Cauchy horizon. The value near the horizon is 
compatible with our previous result, which means that such an approximation may be 
applied to other colliding plane wave space-times. Even with this approximation, the 
initial modes propagated into the interaction region contain a function which cannot be 
calculated exactly and to ensure the correct regularization of the stress-energy tensor 
with the point-splitting technique, this function must be given up to adiabatic order 
four of approximation. 
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1 Introduction 



Although exact gravitational plane waves are very simple time dependent plane symmetric 
solutions to Einstein's equations [jl|, they show some nontrivial global features which are 
essentially due to the non-linearity of General Relativity. One of these features is the absence 
of a global Cauchy surface, which is a consequence of the focusing effect that the waves 
exert on null rays. Another feature, related to the absence of a global Cauchy surface, is 
the presence of a Killing- Cauchy horizon which may be physically understood as the caustic 
produced by the focusing of null rays [0. The inverse of the focusing time is a measure 
of the strength of the wave, thus for an Einstein-Maxwell plane wave such inverse time 
equals the electromagnetic energy per unit surface of the wave. This makes exact plane 
waves very different from their linearized counterparts, which have no focusing points and 
admit a globally hyperbolic space-time structure. One expects that exact plane waves may 
be relevant for the study of the strong time dependent gravitational fields that may be 
produced in the collision of black holes [0, ^] or to represent travelling waves on strongly 
gravitating cosmic strings In recent years these waves have been used in classical general 
relativity to test some conjectures on the stability of Cauchy horizons |^, and in string 
theory to test classical and quantum string behaviour in strong gravitational fields |P, |^, 10 



Their interest also stems from the fact that plane waves are a subclass of exact classical 
solutions to string theory ||lT], |12], 

When plane waves are coupled to quantum fields the effects are rather trivial since they 
produce neither vacuum polarization nor the spontaneous creation of particles. In that sense 
these waves behave very much as electromagnetic or Yang-Mills plane waves in fiat space- 
time [|1^, |I5|. Still the classical focusing of geodesies has a quantum counterpart: when 
quantum particles are present the quantum field stress-energy tensor between scattering 
states is unbounded at the Cauchy horizon, i.e. where classical test particles focus after 
colliding with the plane wave 0]. This suggests that the Cauchy horizon of plane waves 
may be unstable under the presence of quantum particles. The classical instability of the null 
Cauchy horizons of plane waves is manifest when non-linear plane symmetric gravitational 
radiation collides with the background wave, i.e. when two plane waves collide. In this 
case the focusing effect of each wave distorts the causal structure of the space-time near 
the previous null horizons and either a spacelike curvature singularity or a new regular 
Killing-Cauchy horizon is formed. However, it is generally believed that the Killing- Cauchy 
horizons of the colliding plane wave space-times are unstable in the sense that "generic" 
perturbations will transform them into spacelike curvature singularities. In fact, this has 



been proved under general plane symmetric perturbations [|1^. Also exact colliding plane 



wave solutions with classical fields are known that have spacelike curvature singularities 
and which reduce, in the vacuum limit, to colliding plane wave solutions with a regular 
Killing-Cauchy horizon 



Colliding plane wave space-times are some of the simplest dynamical space-times and, as 
such, they have been used as a testing ground for some problems in classical general relativity 
such as the just mentioned stability of the Killing-Cauchy horizons, or the cosmic censorship 
hypothesis [0 . Note that in a colliding plane wave space-time with a Killing-Cauchy horizon 
inequivalent extensions can be made through the horizon, and this implies a breakdown of 
predictability since the geometry beyond the horizon is not uniquely determined by the initial 
data posed by the incoming colliding plane waves. The singularities in colliding wave space- 
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times are also different from the more familiar cosmological and black hole singularities which 
originate from the collapse of matter since they result from the non-linear effects of pure 
gravity. The type of singularities also differs in the sense that these are all-encompassing, 
i.e. all timelike and null geodesies will hit the singularity in the future. 

In three previous papers [|1^, |20|, ^ we studied the interaction of massless scalar quantum 
fields with a gravitational background which represents the head on collision of two linearly 
polarized shock waves followed by trailling gravitational radiation which focus into a Killing- 
Cauchy horizon. The space-time is divided into four distinct regions: a fiat space region 
(which represents the initial fiat region before the waves collide), two single plane wave 
regions (the plane waves before the collision) and the interaction region which is bounded 
by the previous three regions and a regular Killing- Cauchy horizon. Each single plane wave 
region contains a type of topological singularity usually referred as folding singularity which 
is a remnand of the coordinate singularity that a free gravitational plane wave develops as 



a consequence of its focusing properties over null geodesies |22|. The interaction region is 



locally isometric to a region inside the event horizon of a Schwarzschild black hole with the 
Killing-Cauchy horizon corresponding to the event horizon p3|, E3. The presence of the 



Killing horizon made possible the definition of a natural preferred "out" vacuum state [25 



and it was found in ^ that the initial fiat vacuum state contains a spectrum of "out" 
particles. In the long wavelength limit the spectrum is consistent with a thermal spectrum 
at a temperature which is inversely proportional to the focusing time of the plane waves. 
Of course, the definition of such "out" vacuum is not possible when we have a curvature 
singularity (i.e. in the "generic" case) instead of an horizon, whereas a physically meaningful 
"in" vacuum may be defined in all colliding plane wave space-times. 

In reference we computed the expectation value of the stress-energy tensor of the 
quantum field in a region near both the horizon and the folding singularities in the initial 
flat space vacuum. We found, not surprisingly, that the stress-energy tensor is unbounded 
at the horizon. The specific form of this divergence suggests that when the backreaction 
is taken into account the horizon will become a space-time singularity, i.e. the Killing- 
Cauchy horizon is unstable under vacuum polarization. Note that this is a non perturbative 
effect but the result of the nonlinearity of gravity, since gravitational waves in the linear 
approximation do not polarize the vacuum. In fact the vacuum stress-energy tensor of a 



quantum field in a weakly inhomogeneous background was computed by Horowitz p6[, and 



it is easy to see that such a tensor can be written in terms of the linearized Einstein tensor 
only which vanishes for gravitational waves. 

The non perturbative evaluation of the expectation value of the stress-energy tensor of a 
quantum field in a dynamically evolving space-time is generally a difficult task. Even when 
the exact modes of the quantum field equation are known it may not be possible to perform 
the mode sums in order to get the quantum field two-point function or, more precisely, the 
Hadamard function, which is the key ingredient in the evaluation of the stress-energy tensor. 
In our colliding wave space-time we do not even know the exact solution of the modes in 



the interaction region (a similar situation is produced in the Schwarzschild case [^, p9|). 
In the previous paper |2l| we managed to calculate these mode sums using the fortunate 
fact that the geometry of the colliding space-time is such that the initial modes which come 
from the fiat region are strongly blueshifted in their frequency in the interaction region near 
both the horizon and the folding singularities. In the present paper we will deal with the 
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problem of calculating these mode sums for a more generic colliding plane wave space-time. 
Nevertheless, for simplicity, we will consider the particular space-time described above and 
we will discuss the possible generalizations to other colliding plane wave space-times. The 
main goal of the paper will be the introduction of a suitable approximation to calculate 
such mode sums in the causal past of the center of the collision, which will allow us to 
approximately compute the evolution of the stress-energy tensor from the beginning of the 
collision until the formation of the Killing-Cauchy horizon. In fact, in the region close to the 



horizon we will recover the behaviour of the stress-energy tensor computed in ||2T[. We will 



see that, in order to propagate the initial fiat mode throughout the causal past of the center 
of the collision, only a subset of the Cauchy data is required and furthermore the behaviour 
of the differential equations satisfied by these modes may be conveniently approximated. We 
will see that these approximations can be also recovered by a convenient change in the space- 
time geometry throughout the causal past of the center of the collision. In fact, in order to 
evaluate the mode sums involved in the Hadamard function and to correctly substract the 
nonphysical divergences, with the point-splitting prescriptions, it is necessary to work with 
approximations in the space-time geometry rather than in the differential equations. 

The plan of the paper is the following. In section 2 the geometry of the colliding plane 
wave space-time is briefly reviewed. In section 3 the mode solutions of the scalar field 
equation are given for the four different regions of the space-time, it is only in the interaction 
region that exact solutions for these modes cannot be found and therefore we develop a 
suitable approximation in order to find them, first by doing approximations in the differential 
equations and then by conveniently changing the space-time geometry. In section 4 the point 
splitting technique is reviewed and adapted to the computational purposes of this paper. In 
section 5 this technique is used to regularize the Hadamard function and calculate its value 
by a mode sum in the causal past of the center of the collision. Finally, in section 6 the 
expectation value of the stress-energy tensor near the horizon is calculated. A summary and 
some consequences of our results, such as the backreaction problem, the quantum instability 
of the Killing-Cauchy horizon and the generality of these results are discussed in section 7. 
In order to keep the main body of the paper reasonably clear, some of the technical details 
of the calculations have been left to the Appendices. 



2 Description of the geometry 

We will consider in this paper an example of a colliding plane wave space-time where instead 
of a curvature singularity a Killing-Cauchy horizon is produced. It represents a collision of 
two pure gravitational shock waves followed by trailing radiation, with an interaction region 
which is locally isometric to a region of the interior of a Schwarzschild black hole, where the 



Killing-Cauchy horizon corresponds to the black hole event horizon [^, 2^|. The space-time 
contains four space-time regions (see Fig. 1), given by 

ds'jy = ALiL2dudv — dx^ — dy'^, (1) 

1 — Sill 1) 

dsijT = ALiL2{l + ^invYdudv rfx^ — (1 + sinf )^ cos^ f c/y^, (2) 

1 + sin V 

X — sill u 

ds'jr = 4LiL2(l + sin u)'^dudv dx^ — (1 + sin -u)^ cos^ u dy"^, (3) 

1 + sin u 
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ds] = ALiL2\l + sm.{u + v)?'dudv — — r (ix^ — [l + sin(-u + f )1^cos^(m — f) (iu^, (4) 

1 + sm(u + v) 

where for convenience we have used u and v as dimensionless null coordinates, and where Li 
and L2, are length parameters such that L1L2 is directly related to the focusing time of the 
collision, i.e. to the inverse of the strength of the waves, which is a measure of the amount 
of nonlinearity of the gravitational waves [22|. 



This colliding wave space-time, as shown in Fig. 1, consists of two approaching waves, 
regions II and III, in a flat background, region IV, and an interaction region, region I. 
The two waves move in the direction of two null coordinates u and f , and since they have 
translational symmetry along the transversal x-y planes, the interaction region retains a 
two-parameter symmetry group of motions generated by the Killing vectors and dy. The 
four space-time regions are separated by the two null wave fronts u = and f = 0. Namely, 
the boundary between regions I and II is {0 < m < 7r/2, v = 0}, the boundary between 
regions I and III is {m = 0, < v < 7t/2}, and the boundary of regions II and III with 
region IV is {m < 0, v = 0} U {u = 0, v < 0}. Region I meets region IV only at the surface 
u = V = 0. The Killing-Cauchy horizon in the region I corresponds to the hypersurface 
u + V = tt/2 and plane wave regions II and III meet such a Killing-Cauchy horizon only at 
V = {u = 7r/2, V = 0} and V' = {u = 0, v = 7r/2} respectively. Observe that plane wave 
regions II and III contain a singularity at m = 7r/2, for region II, and v = 7r/2, for region 
III. These singularities are not curvature singularities but a type of topological singularity 
commonly referred to as a folding singularity This terminology arises from the fact 

that the whole singularity u = 7r/2 in region II (or = 7r/2 in region III) must be identified 
(i.e. "folded") with V (or V) (see for more details and for a 3-dimensional plot of this 
space-time). 



3 Mode propagation 

In this section we will study the interaction of a quantum field with the colliding plane wave 
background. In particular we will be interested in the value of the quantum field which 
corresponds to the initial vacuum state in the flat region, all over the causal past of the 
collision center (region S in Fig. 2). This turns out to be a geometrical problem which 
consists in solving the field equation in the four space-time regions and smoothly joining the 
distinct solutions. We start with the field solution in the fiat region, which is chosen to be 
the usual vacuum state in Minkowski space-time. This vacuum solution will set a well posed 
initial value problem on the null boundaries S = {m = 0, f < 0}U{m < 0, v = 0}, by means 
of which a unique solution for the field equation can be found throughout the space-time, 
i.e., in the plane wave regions (regions II and III), and in the interaction region (region I). 
However, although it is rather easy to find the solution of the field equation in regions II 
and III which matches smoothly with the boundary conditions, it turns out to be a difficult 
problem for the interaction region. The reason is essentially due to the intrinsic differences 
between the the geometry of the plane wave regions and the interaction region. In fact, the 
plane wave regions are either conformally fiat or type N in the Petrov classification, but the 
interaction region can be more generic. We will refer, from now on, to this problem as the 
mode propagation problem. 

For certain particular examples it has been only partially possible to solve the mode 
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propagation problem |T9|, ^ , essentially due to the separable properties of the field 

equation. There is not, however, an analitical way to solve such a problem in general, for 
any colliding plane wave space-time. Nevertheless, if we are only interested in the collision 
center, we do not need to solve the mode propagation problem throughout the entire space- 
time but only in the causal past of the collision center (i.e region S). We will see in what 
follows that this not only will introduce an important simplification but also will allow us to 
find an approximation procedure, presumably valid for any colliding plane wave space-time. 

Let us now solve the field equation in all four space-time regions. We will consider for 
simplicity a massless scalar field, which satisfies the usual Klein-Gordon equation, 

□0 = 0, (5) 

and let us take the line element, 

which applies globally to the four space-time regions, and where the functions U, V and 
M, can be directly read from (|lD-(^. Then, the field equation can be separated in a plane- 
wave form solution for the transversal coordinates x and y, with fc^. and ky, respectively, 
as separation constants. This plane- wave separation is just a trivial consequence of the 
translational symmetry of the space-time on the planes spanned by the Killing vectors dx 
and dy. The field solution is thus, 

V, X, y) = e^("'")/2 f{u, v) e'^^^+'^yy ^ (7) 
where the function f{u,v) satisfies the following second order differential equation, 

fe-^/2) . 

+ n{u, v)f = 0; Q{u, v) = -^-ZuJT^ + Tf'''^'' {kle-"" + ky) . (8) 

From now on, we will refer to Q{u,v) as the potential term. In the fiat region (region IV) 
this potential term is simply, 

niyiu,v) = LMkl + ^y), (9) 

where we have used that the functions U and V in (^) are zero in the fiat region, and the 
funtion exp(— M) = 2L1L2. Using (|), equation (|) can be solved as, 

= e-'2^+"-'2^-^ (10) 

where k± are two new separation constants with dimensions of energy, and we define for 
convenience two dimensionless constants as k± = y/LiL2 k±. These new separation constants 
are directly related to the previous ones k^ and ky by, 

Ak^k_ = kl + kl. (11) 

The field solution in region IV reduces, thus, to the usual Minkowski plane wave solution, 
i.e, 



5 



(j)k{u,v,x,y) 



2A;_(27r) 



-i2k^u—i2k-V+ikxX+ikyy 



'12) 



These modes are well normalized on the null hypersurface E = {{u = 0, v < 0} U {u < 
0, V = 0}, which is the boundary of the plane wave regions II and III with the flat region 



IV. Even though E is a null hypersurface, a well defined scalar product is given by (see [jT9| 
for details), 



dxdy 



/-O 

du + 

v=0 J —oo 



dv 



u=0 



(13) 



The modes (|T2|) will determine on S a well posed set of boundary conditions for modes in 
regions II and III. There, the potential term in equation (^ is simply, 



-Mi+Ui f u2-Vi 



-e ( kZe + kfe 



(14) 



where the label i = II or i = III in the functions U, V and M, stands for their particular 
values in the plane wave regions II or III. Then, the solution of equation (^ in regions II 
and III with the boundary conditions imposed by the flat modes (|l^) on the hypersurface 
S, can be easily found as. 



-i2k-v-iAu{u)/{2k-) . 



-i2k+u-iAiu{v)/{2k+) . 



in region II, 



in region III, 



where the generic function Ai{(), with i =11, III is given by, 

A,{C) = /^<-e-*^'«')+^'(^') (e-^'«')A;2 + e^»(^')/t2\ 
Jo 2 ^ 

Therefore, the well normalized "in" modes in regions II and III are. 



(15) 



(16) 



^2k_{2Txf 



^ikxX-\-ikyy 



COSM 
1 

( 

COS V 



-i2k- D-i Ajj(n) /(2k- ) . 



in region II, 



-i2k+u—iA,„{v) / {2k+) . 



; in region III, 



where, using the notation ki = ^/LlL^kx:, k2 = y/LiL^ky, 



l + sinC)^(9-sinC) , 15 15 " 

: cosC C — 12 

2cosC 2 ^2^ 



+ kl tan (, 



(17) 



(18) 



Now the Cauchy problem is well posed on the boundaries {u = 0, < v < vr/2} and 
{0 < M < 7r/2, V = 0} between plane wave regions II and III and the interaction region. 
Notice that the initial modes (|T7|) are well normalized on the boundary E between the fiat 
region and the plane wave regions, and this means, from general grounds, that they remain 
well normalized on the boundary between the plane waves and the interaction region. 
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We now have to solve equation in region I with the boundary conditions imposed 
by (|T7|) on the hues Sj = {u = 0, < v < 7r/2} U {0 < m < 7r/2, v = 0}. These 
hues are characteristic hues for the partial differential equation and therefore the only 
independent boundary conditions, i.e. the Cauchy data, are the initial values of the function 
f{u, v) on them (the normal derivatives of the function on the characteristics, which are 
usually part of the Cauchy data, are determined by the values of the function f{u,v) itself 
pT|). In fact, we only need to find the solution of equation (|^) in the neighborhood of the 



collision center and from general grounds in partial differential equation theory (see |3T| for 
details) this solution will be determined only by a subset of the Cauchy data (see Fig. 2). 
We will see in what follows that since that subset of Cauchy data includes only smooth 
non-singular functions, the mode propagation problem can be suitably approximated. We 
start with the change of coordinates, 

t = u + v, z = V — u, (19) 

in equation (§) and we obtain, 

f,u-f,.. + ^{t,z)f = 0, (20) 
where the potential term is given by, 

{l + smtrkl + {smH)/A ^ fcj - (sin^ z)/4 ^ 
' cos^ t cos^ z 

The coordinate singularity (Killing-Cauchy horizon) occurs aX t = u + v = tt/2. The lines 
Si = {m = 0, < f < vr/2} U {0 < m < 7i/2, v = 0}, where the Cauchy data is imposed by 
(P!7D, are characteristics for the equation (|), or (PD|). The collision center is determined by 
the simple condition u = v, and thus the subset of Cauchy data that affects the neighborhood 
of u = V lies on the lines. Si = {u = 0, < v < 7r/4} U {0 < m < 7r/4, v = 0}. Region 
S in Fig. 2 is the causal future of this Cauchy data (or equivalently, the causal past of the 
colision center). 

However, the behaviour of the variables t and z in equation ( PI1[ ) in region S is very 
different. Since t runs from to tt/2 and z runs from — 7r/4 to tt/4 in this region, the 
potential term ( pT]) blows up as coordinate t goes to 7r/2, but is perfectly smooth over the 
entire range of coordinate z. This fact suggests that in the whole region S the physical 
results that we may expect are directly related to the coordinate t and we may not expect 
any physically remarkable change if we take 2; = in equation (|20|) . However, if we want to 
be consistent with such an approximation, we must also modify the boundary conditions that 
lie on the line segments Si. Since on the boundary Si we have that t = ±z and coordinate 
t runs from — 7r/4 to 71/4, we must also take t = z = 0. This means that the boundary 
conditions on Si, given by (p!7|), reduce in such an approximation to the flat boundary 
conditions ([121) . Therefore we change the mode propagation problem for the colliding wave 
space-time into a rather simpler problem, which is clear from Fig. 3, and which requires 
only that we find a solution to equation ( ^0]) with initial conditions given by the Minkowski 
flat modes (|T2p below the hypersurface {t = 0, — 7r/4 < z < vr/4}. In fact, using such a 
simplification we can explicitly eliminate the dependence of equation (|20|) on coordinate z, 
with z = 0, by taking a 2;-plane wave solution, i.e., f(t,z) = g(t)e^^'^^, where is a new 
separation constant. 
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Observe that we have changed a difficult problem of two-variable partial differential 
equations into a rather easy one-dimensional Schrodinger-type problem. Recall, however, 
that none of the discussion above is applicable when a solution of equation (^) in the 
neighborhood of the folding singularities V and V is required. This is not only because 
in that case both coordinates t and z take values near 7r/2 and thus the potential term 



(pTD is unbounded as z 7r/2, but also because the boundary conditions (Jl7\} are also 
unbounded as the folding singularities at t = ±z = 7i/2 are approached. In that case the 
mode propagation problem is much more complicated and a more detailed discussion is 



required (see [|T|, ^ ^ 0) for details). 

In order to solve this Schrodinger-type problem, rather than relying on the discussed 
approximations for the exact field equation (|20|) , we will rewrite a new field equation using 
an adequate approximation for the line element throughout the entire causal past of the 
collision center. This new approach, which may seem redundant in the case of the particle 
production problem, is absolutely necessary when the renormalization of the stress-energy 
tensor is discussed. This is essentially because the process of renormalization involves the 
subtraction of the infinite divergences that arise from the formal definition of the stress- 
energy tensor, and these divergences can be expressed as entirely geometric terms, which 
are independent of any possible approximations in the field equation. This means that in 
order to recover the geometric divergences in the stress-energy tensor, any approximation in 
the field equation solutions must be related to a suitable approximation in the space-time 
geometry. 

Approximating the field equation (|20|) in the causal past of the collision center by taking 
z = is essentially equivalent to changing the line element, in the causal past of the collision 
center, by a related line element obtained directly from the original (|^) by setting z = 0, 
i.e., 

dsl = LiL2(l + sint)^ (rft^ - dz^^ - j ~ ^^^^ - (1 + smtfdy"^. (22) 



We will suppose that the line element (|22D applies all over the causal past of the collision 
center, not only in the interaction region but also through the plane wave regions II and III 
in the sense of Fig. 3. The plane wave collision starts at t = but to avoid smoothness 
problems derived from such an approximation, we will suppose that (|2^) applies exactly on 
a range e < t < n/2, for a certain e > 0. In the range < t < e, as described below, 
we will interpolate a line element which smoothly matches with the fiat space at t = 0. 
Nevertheless, the details of this matching will not affect the main physical features. 
The exact field equation for this approximate space-time is, 

(□ + ei?)0 = O, (23) 

where it is necessary to consider a coupling curvature term in the field equation because, 
although the exact space-time is a vacuum solution, we have a bounded nonzero value for R 
in the approximated space-time. In order to solve this new field equation, we start rewriting 
the line element (^) in the following general way, 

ds' = ihhh) df^ - dz- _ j^:^^ dx^ _ dy\ (24) 
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where the are functions of coordinate t alone, which for values of0<e<t<7r/2, can be 
straightforwardly determined by direct comparison with (^) as fi{t) = \fL^\L^(\ + sint)^, 
/2(t) = y/LiL^ cost, fs(t) = cost. For values t < we take fi(t) = f2(t) = ^/L^I^, 
f3{t) = 1, which correspond to their values in flat space. Finally, in the interval < t < e, 
we smoothly interpolate each fi(t) {i = 1,2,3) between these values. Also, in order to 
prevent singularities in the field equation, we conveniently reparametrize coordinate t, by 
t*(t), as follows, 

dt* _ 1 

l^^WY ^^^^ 
Now, we use the following ansatz for the field solutions, 

0fc = /i(r)e^^^^+^'^^^+^'=^", (26) 

where the plane wave factor in coordinates x, y is related to the translational symmetry of 
the space-time along the transversal directions x, y, and the plane wave factor in coordinate 
z is just a consequence of our approximation. Then equation (|23| ) directly leads to the 
following Schrodinger-like differential equation for the function h{t*), 

h,n* + u:\t) h = 0, V{t) ^ u:\t) = f^{t) + fl{t) kl + /|(t) fcj + fl{t) kl (27) 
where the function /o(t) stands for, 

m = [/i(t)/2(t)/3(t)] ^R. (28) 

Such differential equation can be WKB solved, essentially because the short wavelength con- 
dition holds, i.e. u^^d/dt* Ina; <C 1. Observe that this condition reduces to (dt/dt*) dV/dt -C 
2uj^, which becomes particularly accurate when the Killing-Cauchy horizon is approached 
since in that case dt/dt* = fsit) 0. Therefore, the mode solutions (pk which reduce to the 
fiat mode solutions in the region prior to the arrival of the waves, are 

^ _^ikxX-\-ikyy-\-ikzz-i J W{C)dC (29) 



^{2Tr)^2k_W{t) 

where we denote u"^ = kf + k^ + kl with ki = a/ -^1-^2 kx, /c2 = \/L\L2 ky, k^ = kz and 
where W{t) stands for an adiabatic series in powers of the time- dependent frequency u(t) 
of the modes and its derivatives. Up to adiabatic order four (i.e. up to terms involving four 
derivatives of 00 (t)) W{t) it is given by, 

^^,/ ^ ^2 Bo Aa Ba Ca Da , , 

W{t) = ^ + ^ + ^ + 4 + 4 + 4 + -i, (30) 



00'^ UJ^ UJ^ Lj'^ UJ^ UJ^^ ■ 



where, using the notation V = dV/dt*, 

^ 32' ^ 128 ' 258 ' 2048 ' 
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and An, Bn, ... denote the n adiabatic terms in W(t). Up to adiabatic order zero it is 
simply W{t) = uj{t). Observe that in the fiat region prior to the arrival of the waves we 
have W{t) = Cj = (/c^ + + Observe also that since /s = 1 in this fiat region, we 

can use (^5]) to set t* = t, where without loss of generality we choose the origin = at 
t = 0. Therefore, the mode solutions in the fiat region reduce to. 



iIV ~ ikxX+ikyy+ikzZ—iujt 



(27r)32A;_ 

which indeed are the fiat mode solutions defined in (p!2D, recalling that the new separation 
constant kz = is related to the original k± by the ordinary null momentum relations, i.e., 

Cj = kj^ + k_, kz = k+ — k_. (32) 

It is important to understand that we are constructing a set of mode solutions as an 
adiabatic series in terms of derivatives of the frequency uj{t) in the differential equation (p^). 
This procecure is similar but not equivalent to the construction of an adiabatic vacuum state 
where the field modes are expanded as an adiabatic series in terms of the derivatives of 
the metric coefficients (see for example ||3^ for details). In fact, observe for instance that 



the term /g (^) in ( JI7\ j involves two derivatives of the metric since it is directly related to 
the curvature scalar and therefore it would be an adiabatic order two term for an eventual 
adiabatic vacuum construction, but it is simply an order zero term in our adiabatic series 
in derivatives of Lj{t). 



4 "Point splitting" regularization technique 

In this section we briefly review, for the computational purposes of the following sections, the 
"point-splitting" regularization technique to calculate the expectation value of the stress- 
energy tensor of a scalar quantum field in some physical state. The stress-energy tensor 
of the field may be obtained by functional derivation of the action for the scalar field with 



respect to the metric. When the field is massless, it is [34 



T^, = (1- 20 + (2^-^) (7m. (7"^</);«</);;3- 2^0;^. + ^e^7/..</>V"V,^ 



-e i^R^u - -Rg^^u + i^^Rg^uj (t>\ (33) 

where ^ is the coupling parameter of the field to the curvature. 

To quantize, the field is promoted into a field operator acting over a given Hilbert 
space [|^, |3^, 0(x) = Y.k a-kUkix) + a\u%{x), where the standard creation and 



annihilation operators and {uk{x)} is a complete and orthornormal set of solutions of the 
Klein-Gordon equation (^3)). Mathematically the field operator (f){x) is a point distribution, 
therefore, the quantum version of the stress-energy tensor ( PB| ) is mathematically patholog- 
ical because it is quadratic in the field and its derivatives. One possible way to give sense 
to that expression is to note that the formula (^) can be formally recovered as. 
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(T^,(x))= \imV^,G^^\x,x'), (34) 

where x') is a Green's function of the field equation defined as the vacuum expectation 

value of the anticommutator of the field, and called the Hadamard function, 

= {{ct>{x),<P{x')]) = Y.{uk{x)ul{x') + Uk{x')ul{x)] . (35) 

k 

As a product of distributions at different points this is mathematically well defined. The 
differential operator V^y is given in our case by, 



V 



:i - 2 ^ (V^. V, + V,. V^) + (2 e - ^) 9^.. \ (V„' V" + V« V°' 



(V^V, + V^/V,0 + g^y I (V„V" + V„/ V°' 



^ R 



(36) 



However, the above differential operation and its limit have no immediate covariant meaning 
because G'^^^x^x') is not an ordinary function but a hiscalar the differential operator 
is nonlocal] thus we need to deal with the nonlocal formalism of bitensors (see, for 



example [^, or the Appendix B of reference ||2T]] for a review on this subject). 

The above procedure still leads to a divergent quantity since we know that even in flat 
space-time G^^\x,x') has a short-distance singularity and that a "vacuum" subtraction has 
to be performed to G^^\x, x') in order to obtain a regularized value. To regularize we assume 
that {x,x'), has a short-distance singular structure given by 



S{x, x') 



+ f (x, x') In (t(x, x') + w{x, x') 



(37) 



(47r)2 L a{x,x') 

where a{x,x') = (1/2) s'^{x,x') is the geodetic biscalar (being s{x,x') the proper distance 
between x and x' along a non-null geodesic connecting them), A is the Van Vleck-Morette 



determinant ||3^, which is singularity free in the coincidence limit, and where v{x,x') and 



w{x, x') are biscalars with a well-deflned coincidence limit for which we assume the following 
covariant expansions. 



t;(x, a;') = ^ f;(x, x')cr'(x, x'), w{x, x') = '^wi{x, x')a''{x, x'). (38) 

1=0 1=0 

A Green's function expressed in this form is usually called an elementary Hadamard solution, 
the name of which comes from the work of Hadamard on the singular structure for elliptic and 
hyperbolic second order differential equations. Note, however, that this Hadamard singular 
structure is not a general feature of any Green's function of the Klein-Gordon equation. In 
other words, although for an extensive range of space-time and vacuum states, the vacuum 
expectation value of the anticommutator of the fleld, G^^\x,x'), has this singular form, 
this is not a general property. However a theorem states that if G^^\x,x') has the singular 
structure of an elementary Hadamard solution in a neighbourhood of a Cauchy surface of 
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an arbitrary hyperbolic space-time, then it has this structure everywhere p5| , |35| . As a 
corollary of this theorem, G^^^ {x, x') has this singular structure if the space-time is fiat to 
the past of a space-time Cauchy surface, as is the case of our colliding plane wave space-time. 



This and other considerations led to a proposal by Wald |^6| that any physically reasonable 



quantum state must be a Hadamard state, that is to say, a state for which G'^^\x,x') takes 
the short-distance singular structure of an elementary Hadamard solution. 

The coefficients Vi{x,x') and Wi{x,x') can be directly obtained by substitution in the 
differential equation, {O^ + ^R)S{x,x') = 0. Recursion relations for Vi{x,x') and Wi(x,x') 
are then obtained These relations uniquely determine all the Vi{x,x') coefficients 

but the coefficients Wi{x,x') can be written in terms of an arbitrary term wq{x,x'). Up to 
order a, v{x,x') is given by 

f (x, x') = ai{x, x') — - a2{x, x') a + ■ ■ ■ , (39) 



where ai{x,x') and a2{x,x') are the Schwinger-DeWitt coefficients plf whose midpoint ex- 



pansion are given, up to the order required for the regularization calculations, by 

ai{x,x') =af^ + af\^a''a^ + ■■■, 02(0;, x') = 4°^ + ■ ■ ■ , (40) 



and the terms and 02°^ are written in Apendix A in terms of geometrical quantities. 



Only the coefficients Vi{x,x') are related to the singular structure of G^^\x,x') in the 
coincidence limit, and they are uniquely determined by the space-time geometry. This 
means that given any two Hadamard elementary solutions in a certain space-time geome- 
try, both have the same singularity structure in the coincidence limit; therefore given two 
vacuum Hadamard states, |0) and |0), G^'^\x,x') = (O|{0(x), 0(x')}|O) and G^ \x,x') = 
{0\{(l){x),4>{x')}\0), they have the same singular structure. Their finite parts, however, may 
differ because the two vacuum states are related to different boundary conditions, which are 
global space-time features. Mathematically this comes from the fact that the term wo{x,x') 
in the elementary Hadamard solution is totally arbitrary; fixing wq{x,x') we fix a particular 
boundary condition. This suggests a possible renormalization procedure |^ |39|, we can 
eliminate the non-physical divergences of any {x, x') without alterations in the particu- 
lar physical boundary conditions by subtracting an elementary Hadamard solution with the 
particular value wo{x, x') =0, which is the value that corresponds to the flat space case. In 
other words, we define the following regularized biscalar, 

G%\x,x') = G^'\x,x') - 5(a;,x')L„=o- (41) 

This particular value for the elementary Hadamard solution, which we may refer to as the 
locally constructed Hadamard function as opposed to its non-local counterpart G^^\ can be 



easily calculated. In fact, setting wo{x,x') = 0, the biscalar w{x,x') reduces to 



w(a;,x') = -^af a + (42) 
and using the following midpoint expansion for the Van Vleck-Morette determinant we find, 

A'/\x, x') = l + A(2)^, a'^a' + A^^^^p^-^ a'^a^^a^a^ (43) 
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where the coefficients A^^-'^p, A^^^ pppf are written in Appendix A. Then the midpoint ex- 
pansion for the locally constructed Hadamard function ( P7D reduces to, 



a 



ln(/i ^(t) 



1 (0) 

-On 'a 

2 2 



(44) 



where we have included an arbitrary length parameter /i in the logarithmic term which, 
as we will mention below, is related to the two-parameter ambiguity of the point-splitting 
regularization scheme. 

Then by means of Gg^ {x,x') we can construct a (^/f/) by differentiation with respect to 
the nonlocal operator (0). This regularization procedure, however, fails to give a covariantly 
conserved stress-energy tensor essentially because the locally constructed Hadamard function 
(^^ is not in general symmetric on the endpoints x and x' (i.e. it satisfies the field equation 
at the point x but fails to satisfy it at x'). In fact, it can be seen that for a massless 
conformal scalar field (i.e. ^ = 1/6), 



1 



1 



V^(Tfj = - lim V. a^, + - R(x') G^(x,x') -- 



1 



V^af (x) 



2 A* 



(45) 



where a^\x) is the coincidence limit of the Schwinger-DeWitt coefficient a2{x,x') given 
in Appendix A. Thus to ensure covariant conservation, we must introduce an additional 
prescription: 



X) 



(46) 



Note that this last term is responsible for the trace anomaly in the conformal coupling 
case, because even though (T^(a;)) has null trace when ^ = 1/6, the trace of {T^^{x)) is 

given by (T^) = —a2'\x)/{lQn'^). Recall that for scalar fields with non-null mass we may 
regularize the Hadamard function by subtracting a truncated DeWitt-Schwinger series (see 
p4| for details). This gives a symmetric regularized Hadamard function by means of which 
an automatically conserved stress-energy tensor can be obtained. Nevertheless, the DeWitt- 
Schwinger series is not well defined in the massless limit and therefore it is not suited to our 
case. 

The regularization prescription just given in ( ^^ satisfies the well known four Wald's 
axioms pBI, HOl, ETl E2|, a set of properties that any physically reasonable expectation value 



of the stress-energy tensor of a quantum field should satisfy. There is still an ambiguity in 
this prescription since two independent conserved local curvature terms, which are quadratic 
in the curvature, can be added to this stress-energy tensor. In particular, the ^u-parameter 

Such a two-parameter 



ambiguity in (P^ is a consequence of this (see for details), 
ambiguity, however, cannot be resolved within the limits of the semiclassical theory, it may 
be resolved in a complete quantum theory of gravity |^ . Note, however, that in some sense 
this ambiguity does not affect the knowledge of the matter distribution because a tensor of 
this kind belongs properly to the left hand side of Einstein equations, i.e. to the geometry 
rather than to the matter distribution. 
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5 Hadamard function in the interaction region 



Here we calculate the Hadamard function G^^'^ (x, x') in the interaction region for the initial 
vacuum state defined by the modes 0^, (pQl). The Hadamard function can be written as, 



G^^\x,x') = Y.<Pk{x)<Pl{x') +C.C. (47) 

k 

Note that solutions 0^ contain the function h(t*), which cannot be calculated analytically 
but may be approximated up to any adiabatic order as described in (|30|)-(|3lD. This means 



that we have the inherent ambiguity of where to cut the adiabatic series. However, observe 
from (|25|) and (|28|) that since dt/dt* and V{t) = uj'^{t) —>■ 16 kf towards the horizon, 
the adiabatic series ( pOf ) reduces to ~ a; near the horizon. This means that we could cut 
the adiabatic series (^) at order zero if we were interested in a calculation near the horizon. 
However, this is only partially true. In fact, it would be true if we were only interested 
in the particle production problem (see |]19|) for details) but it is not sufficient for the 
calculation of the vacuum expectation value of the stress-energy tensor. This is because G^^^ 
calculated with h{t*) at order zero does not reproduce the short-distance singular structure 
of a Hadamard elementary solution (^7|) in the coincidence limit x — ^ x'. The smallest 
adiabatic order for the function h(t*) which we need to recover the singular structure of G^^^ 
is order four, basically because our adiabatic construction af the mode solutions is similar 



(but not equivalent) to an adiabatic vacuum state (see p3| for details). 

Although expanding the function h{t*) in (|47|) up to adiabatic order four will give an 
accurate value for the stress-energy tensor near the horizon, it will also give a suitable 
approximate value for this tensor all over the causal past of the collision center (region S in 
Fig. 3). The reason is that even though the short- wavelength condition, i.e. uj~^d/dt* Inu; -C 
1, is particularly accurate near the horizon it also holds throughout region S. 

In the mode sum ( ^7|) we use the shortened notation Y.k = Io° ^^-/^- S^oo^^x jT'oo^K 
or equivalently J2k = (-^1-^2)"^ c^^i /.^ (^^2 /^'ifes/tt', where the change of variables 
(|32D and the usual notation ki = 1/L1L2 kx, /c2 = V-^1-^2 ky have been used. Therefore we 
have, 



I roo roo roc dki dk2 dk^ 



GW(^X,x') = / / / ""2 »"-3 ^ 

2(27r)3LiL2 J-^J-^J-^ Jw{t)W{t') 



e 



-iJ^^,Wi()dC+ika:ix-x')+iky(y-y')+ik^{z-z') _|_ ^ ^ 



Let us start now with the point splitting procedure. We assume that the points x and x' 
are connected by a non-null geodesic in such a way that they are at the same proper distance 
e from a third midpoint x. We parametrize the geodesic by its proper distance r and with 
abuse of notation we denote the end points by x and x', which should not be confused with 
the third component of {t, z, x, y). Then we expand, 

= x^'{t) = x'^(f + e) = a;^ + e + - a;^ + ■ ■ ■ , (49) 

x'i" = x^(r') =x''if-e) =x^-ex1 + -xi^ , (50) 
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where we have defined, 



Xq 



/I 



dx^ 
dr 



x'i 



d'^x'^ 
~d^ 



(51) 



i.e. the subscript of a given coordinate indicates the number of derivatives with respect to 
r at the point f. The geodetic biscalar can then be written as 



a 



= 2ex^, x^xi^ = E, (52) 

with E = +1 if Xj* is timehke and E = — 1 if Xj* is spacehke and where is the geodesic 
tangent vector at the midpoint x with modulus the proper distance on the geodesic between 
X and x' (see for instance ||38| , 45 1 or the Appendix B of for a review on bitensors). 



Using dH) and (|0D we can write !l,W{C)dC and [W{t)W{t')\~^''^ , in (|8D, in powers 



of e. 

a _ 



Note also that 



X 



2ex\ + (e^/3)x3 + (e^/60)x5 + 0(e'^), where we denote 



X = X, X = y, X 
in powers of e as 



z. With these expressions we can expand the exponential term in (^ 



^-i J^,W{0dC+ikx{x-x')+iky{y-y')+ik3{z-z') _ ^-ie{So uj-Si k:^-52 ky-Sg kz) ^ _|_ . . 

where we do not expand the zero adiabatic terms because they will be useful as new inte- 
gration variables. These terms have been included in the coefficients 6q and 6i as, 

2 4 2 4 

Sq = 2tl + -t; + —t* 5i = 2x\ + -xi + —xi. (53) 
1 3 3 gg 5' * 1 3 3 gg 5 V ; 

Since the number of derivatives d/dt* determine the adiabatic order, we introduce a new 
parameter T which will indicate the adiabatic order. Then at the end of the calculation we 
will take T = 1. With abuse of notation, we also denote, u = uij), W = W{t). Recall that 
in (^) W is given up to adiabatic order four, thus we have now, 

W = uj + T [ — + — ]+T — + — + — + — , 

where the coefficients An, Bn, Cn, are given in (P7| ) and (^Tj). 

Now, if we substitute these expansions in ( ^8]) and separate the different adiabatic terms 
by expanding in powers of the parameter T, we can easily see that G^^\x, x') can be written 
as a sum of the following type of integrals 



OO /"OO 



dk^ dky dk^ 



^-ie(dooJ—di kx—d2 ky~ds kz) (54) 

OO./— OOv'— OO UJ"" 

where, V{t) = fo{t)'^ + kl fi{t)'^ + k^ f2{t)^ + k1 f^it)'^. For computational purposes, these 
integrals may be solved as follows. First we use the following change of coordinates, 

zi = fi k^, Z2 = /2 ky, Z3 = /s k^, dk^dkydkz = {fif2fz)~^dzidz2dz3, 
with fi = fi(t) and the following definitions, 

Ao = So, \ = SJ-\ = fl + zl + zl + zl 
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We can then use spherical coordinates, zi = rsin^cos0, Z2 = rsin^sin0, = rcosO, 
and the following simplifying notation Qi = sin ^ cos 0, ^2 = sin ^ sin 0, ^3 = cos^. Thus, 
= /q + and we conveniently define. 



Observe that we can also eliminate the factor r^, which appears in after the change to 
spherical coordinates, using the differential operator {ie)^^d/dA. We may finally turn the 
terms V^^\t), which appear in into the differential operators. 



[hitff 



-2 



After all these changes, the integral (^J) transforms into 



-1 



dd sin 6 



9A2 



/3 9 



9A2^"' 



where for simplicity we have conveniently defined. 



T 



°° dr 
— I 



-ie(Ao n-Ar) 



(55) 



and the solution of these X„ integrals are discussed in Appendix B. 

Solving integrals (^) and then using (0), the relation between the geodesic co- 

efficients and in terms of t* and x*^ (see Appendix C) and after performing a rather 
tedious calculation, the Hadamard function reduces to. 



G«(x,x') 



112 _ a'^a"" 



-1-2A(2) 



-2A(^) 



fMUpT 



a 



a 



fJ,U 



(56) 



where L is a logarithmic term defined as L = 27 + ln(cr ^i?/2), being 7 Euler's constant and 
where all the involved coefficients A, b, C^^... are given in Appendix A. 

According to (0), the Hadamard function can be regularized using the elementary 
Hadamard solution (^) and finally the regularized expression for G^^\x,x') up to order 



IS, 



-af L 



(af)A( 



2) 



fiu "1 flu 



a^a On cr 

2 ^ 



(57) 



where L is a bounded logarithmic term given by, L = 27 + ln(/i^ ^R/2), ^ being the arbitrary 
length parameter introduced in (0), and where the coefficient B is B = b + 3 a2'^''/(327r^), 
which is given also in Appendix A. From (|57D we can directly read off the regularized mean 
square field in the "in" vacuum state as (0^) = A/2 — af'^L/ (167r^). It is important to remark 
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that the term -D^^^j in (|57|) appears only as a consequence of our approximate procedure 
of calculating the Hadamard function, i.e. using an adiabatic order four expansion for the 
initial modes in powers of the mode frequency uj{t) and its derivatives. Had we used an exact 
expression for the initial modes (or an adiabatic vacuum state |3^), such a term would not 
appear. 

6 Expectation value of the stress-energy tensor 

To calculate the vacuum expectation value of the stress-energy tensor we have to apply the 
differential operator (|36D to (^). As we have already pointed out, this is not straightforward 
because we work with nonlocal quantities. Note first that the operator ( PB]) acts on bitensors 
which depend on the end points x and x', but the expression ( |57D for G^"* depends on the 
midpoint x. This means that we need to covariantly expand (^) in terms of the endpoints 
X and x'. Consider the general expression, 

{x,x') =A + aB + C^^ a"cT^ + D^^^j + L E^-^ a^a'^, (58) 

where it is understood that A, S, C^^, -Dq/j^j, E^^p are functions that depend on the endpoints 
X, x' but are evaluated at the midpoint x. Then, following the formalism of Appendix B of 
reference [^, we can expand in the neighbourhood of x as. 



where before expanding, a homogeneization of the indices by the parallel transport bivector 
^p, i.e, = o"^, has been applied. The covariant expansions at x are given by, 

A = A-^A,a" + ^A,^a°a^ + 0(e=^), B = B + 0{e), (60) 

Now we can apply the differential operator (|36D to (|59D and we will consider a confor- 
mal coupling (^ = 1/6), which, from general grounds, should provide a good qualitative 
description for photons. If we introduce the operators, 

4'-^ = V. + V,. v^, = v^.v,. + v,v^, (61) 

the operator (|36|), for the conformal case, is 



-T)«=l/6 ^ J_ 



By the properties (^) of the geodetic interval bivector a'^ we can prove the following iden- 
tities, 



hm, a = - hm, C^^} a = -2 g,,, lim, aV^ = - lim, aV^ = -4 6^^ 6^^, 
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lim£« = -lim42j 

a;— >a;' ^ (j x^x' ^ (j 



where (■ ■ ■) is the usual symmetrization operator. Note that by a straightforward appUcation 



of Synge's theorem (see for example references [|38|, ^] or Appendix B of reference ||2T|), the 
coincidence limits of C^^^ and differ in a sign when they are applied on bitensors for which 
the first covariant derivative has null coincidence limit. But when they are applied on the 
biscalar A in ( |60D they coincide: lim^^x' ^^i^} A = lim^^x' ^^f}} A = The application 

of the opperators ( ]6T| ) over quartic products of a'^ gives path dependent terms when the limit 
X ^ x' is performed. Therefore an appropriate averaging procedure is required. The most 
elementary averaging is called four dimensional hyperspherical averaging and consists 
in giving the same weight to all the geodesies which emanate from x as follows. First one 
analytically continues to a Euclidean metric the components of the tangent vectors to the 
geodesies which emanate from x. Second, one averages over a 4-sphere, and third, the results 
are continued back to the original metric. It is not very complicated to find the following 
averaging formulae. 



(o-^o-z.) = 3t^u, {<y i^cr uCTrCr p) = -cr'^g^^(ugTp), {a ^a^arCr pa^ar^) = -a^gp(^^gr(pg^n))- 

Note that from the symmetry of the averaging procedure, the average of an odd number of 
components vanishes identically. 

We now have all we need to evaluate the vacuum expectation value of the stress-energy 
tensor by application of the differential operator (|36| ) to the expression ( |59D of G^^\x,x'). 
Observe for example that P^^^/^a = 0. Then, in the orthonormal basis 6i = gH"^ dt, 62 = 
gll'^ dz, 6*3 = gl!^ dx, 9^ = g^!^ dy^ using the trace anomaly prescription (^Bj), we obtain the 
following expectation values (Tpv) in the conformal coupling case and for values < e < t < 
7r/2 of coordinate t 



27 + ln(/i^e^/2) ^. 1 
(T,.) = 2880 {L^L^f (1 + z)^ ^"g^"^' " 4 ^^^^ + 

diag (pi it) , -p2 it) , pi it) + 2p2 it) , -p2 it) ) , (62) 

where we have used for simplicity the notation z = sint. The trace anomaly term in that 
case is given by (|4^) as, 

l{K}g,u --^9,u- 2880(LiL2)2 7r2(l + ^)^ ^""^^ ' " ' ' ' 
and the functions pi(t) and P2(^) are given by, 

_ 30254 - 106858 z + 168510 z^ - 147520 z^ + 62495 z^ 

~ 241920 (LiL2)2 7r2 (1 - ^)2(1 + 2)6 ' ^^^^ 



18 



_ -566190 + 1821200 z - 2415234 + 1710444 z^ - 472795 
^'^^^ ~ 2661120 {UUf 7r2 (1 - 2)2(1 + zf ' ^ ^ 

Recall that is a positive definite function in the interval 0<e<t<'7r/2. Both 

functions are unbounded at the horizon it = 7r/2), and the expectation value of the stress- 
energy tensor near the horizon is approximately given by, 

(r^.)|,,./2 = diag(p(t), -Pit), 3p{t), -Pit)) - l{T;)g,,, (65) 

where p{t) = A(LiL2)^^ cost^^ and A ~ 0.00018. This result is compatible with our previous 
result ||21|. In that case, however, we found a bigger value for the parameter A ~ 0.001. 
Such a difference makes sense if we recall that in we obtained an approximate value 
for the expectation value of the stress-energy tensor, essentially using the "blueshift" effect 
as the initial modes approach the folding singularities of the interaction region. Therefore, 
such a result was a good approximation in the region near both the folding singularities and 
the Killing-Cauchy horizon, where a priori we would expect a stronger effect than in the 
center of the collision, which is the only region considered in the present paper. 

For values of t < 0, (T^u) = 0, and to be consistent with the approximation we have 
used for the space-time geometry, we should require that the value of (Tf^u) (§3)' "which is 
valid for0<e<t<7r/2, goes smoothly to zero as t — 0. In fact, this can be achieved 
using an adequate matching of the line element (^) with the flat line element through the 
interval < t < e. Observe that the coefficient of the logarithmic term in the stress-energy 
tensor ( |62t ) deppend only on locally constructed curvature terms (as can be seen from the 
Hadamard function (|57|)). Therefore, with an adequate matching of the space-time geometry, 
this coefficient will also smoothly vanish towards the fiat space region, below t = 0. The 
details of such a matching, however, will not affect the main features of the stress-energy 
tensor (|6^) , particularly when the Killing-Cauchy horizon is approached. 

We must recall, however, that although the value (p2| ) for {Tf^u) satisfies asymptotically 
the conservation equation close to the Killing-Cauchy horizon, it does not satisfy exactly the 
conservation equation throughout region S, essentially because it is obtained by means of an 
approximation in the field modes. Nevertheless, we could obtain a truly conserved (Tfj^u), in 
the context of the present approximation, by solving the conservation equation considering 
a (T^^) given by (T/(t), T|(t), T^^(t), T^^(t)) with the conditions: i) Ty{t) = T^{t), which 
not only is compatible with ( |62D but also holds for the exact interaction region geometry, ii) 
trace anomaly condition, i.e. T^(t) = (T^) — Tl{t) — 2T|(t), iii) the ansatz T/(t) = pi(t), 
which is the approximate value of T/(t) obtained in our calculation. Finally, the conservation 
equation gives straightforwardly a value for T^{t) and T^(t) which are compatible with the 
value Tj^(t) = p2(t) andr^'(t) = —pi{t)—2p2(t) obtained in our approximation. In particular, 
they have the same behaviour near the Killing-Cauchy horizon. 



Inspection of (^) shows that not only is the weak energy condition satisfied |^ , which 
means that the energy density is nonnegative for any observer, but also the strong energy 
condition is satisfied. 
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7 Conclusions 



We have calculated the expectation value of the stress-energy tensor of a massless scalar field 
in a space-time representing the head on collision of two gravitational plane waves through- 
out the causal past of the collision center and in the field state which corresponds to the 
physical vacuum state before the collision takes place. We have performed the calculations 
in this particular region essentially because we could introduce a suitable approximation 
to the space-time metric (see Fig. 3) which not only allowed us to dramatically simplify 
the calculations but also to keep unchanged the main physical features, in particular the 
behaviour of the stress-energy tensor near the Killing- Cauchy horizon of the interaction re- 
gion. In fact, such an approximation is also valid for more generic plane wave space-times, 
and this will be the subject of a forthcoming paper. 

The results we have obtained are the following: before the collision of the waves {T^J) = 0, 
which correspond to the lower edge Fig. 3. Then, after the collision the value of (T^j^u) starts 
to increase until it grows unbounded towards the Killing-Cauchy horizon of the interaction 
region. The weak energy condition is satisfied, the rest energy density is positive and diverges 
as cos~^ t. Two of the principal pressures are negative and of the same order of magnitude of 
the energy density. The strong energy condition is also satisfied, (T^) is finite but {T^^'^) {T^^) 
diverges at the horizon and we may use ref on the stability of Cauchy horizons to argue 
that the horizon will aquire by backreaction a curvature singularity. Thus, contrary to 



simple plane waves, which do not polarize the vacuum []I4|, |15|, the nonlinear collision of 
these waves polarize the vacuum and the focusing effect that the waves exert produce at the 
focusing points an unbounded positive energy density. Therefore, when the colliding waves 
produce a Killing-Cauchy horizon, that horizon is unstable by vacuum polarization. 

In the more generic case when the wave collision produces a spacelike singularity it seems 
clear that the vacuum expectation value of the stress-energy tensor will also grow unbounded 
near the singularity. In fact, in a forthcoming paper we will extend the approximation 
introduced in the present work to a more generic plane wave spacetime with the objective 
of more generally proving that the negative pressures associated to the quantum fields could 
not prevent the formation of the singularity. 
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A Useful adiabatic expansions 

The coefficients for the Hadamard function in (^), using for simplicity the notation z = sin t, 
are: 
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A 



B 



Q + z 
60LiL2 7r2(l + z)3' 

1859 - 9652 z + 17478 z"^ - 14644 z^ + 5543 z^ 
60480 (LiL2)2 7r2 (1 - ^)2(l + zf 

-20959 + 1075 z + 122490 z"^ - 174964 z^ + 76711 
166320 L1L2 7r2 (1 - 2)2(1 + 

30629 - 127830 z + 229803 _ 212506 z^ + 88518 z^ 
332640 (LiL2)2 7r2 (1 - ^)(1 + zf 



D 



yyyy 



D 



zzyy 



L1L2 

282 - 1039 z + 1656 z^ - 1407 z^ + 516 z^ 
3840 7r2(l -2)2(1 + ^)2 

1 + 10 z - 14 ^2 _ 6 ^3 + 39 ^4 

960 (LiL2)2 7r2 (1 + 2)8 



1 



(LiL2)2 ^^^^^ 

-83 + 108 2 + 190 ^2 - 426 z^ + 237 
5760LiL2 7r2(l-z)(l + 2)5 

-184 - 183 z + 1780 2^ - 2359 z^ + 954 z^ 
11520 LiL2 7r2(l-z)2(l + 2)2 



xxyy 



L1L2 



The coefficients for the midpoint expansion of the locally constructed Hadamard function 



(H) are: 



.(0) 



,(2) 



A(2) 

1 /I 



180 



0/3 



-f- 

24 VIO 



pa/3 p _|_ po/37 p 

3 1 1 ^ 



180 



p 



a/375 ) 
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B Some useful integrals. 



Let us define 



/oo 
dk e-''^^° ^\ n = ypTA; A > 0. 
-oo 



(66) 



For the particular value n = —1 this integral can be easily solved. Performing the following 
change of variable, 



t = t{k) = Ao f2 - Ai = Ao VWTX - Ai k, 



(67) 



we must consider separately the two possibilities |Ao| > |Ai| or |Ao| < |Ai|. For |Ao| > |Ai| 
and without loss of generality we may take Aq, Ai > 0, and invert the change (|67|) as, 



k = k{t) = 7 



where we use two new variables 7 and /?, defined as. 



7 



-1 



A^ 



{6t 



(69) 



Note that k{t), in (p5|), is a bivalued function and we have to be careful in changing the 
integration limits. Since |Ao| > |Ai|, then \imk^±oot{k) = +00, and this means that ( |67D 
has an absolute minimum, which is at the point (+7/3Ai, (3). With this it is easy to see 
that we have to take as inverse function of ( |67D the function (|^) with the plus sign to the 
right of the minimum t = j3 and with the minus sign to the left of the minimum. Therefore 
we can split the integral (B§), for n = —1, into two parts on each side of t = f3, i.e. 



7/3Ai dk 
-00 



g-ie(Ao fJ-Ai fc) _|_ 



°° dk 



-ie{Ao n-Ai k) 



(70) 



/7/3A1 Vt 

The change of variables ( p7D gives Vt'^dk = ±(t^ — P'^)~^^'^dt, with the minus sign for the 
first integral in ( ffOD and the plus sign for the second integral. Then the integration can be 
easily performed in terms of a zeroth order Bessel function El], as 



-ie t 



■.dt = 2Ko {ie (3) . 



(71) 



For the case |Ao| < |Ai|, we define the parameter 7, as, 7^^ = A^ — Ag, and the inverse 



function of (E7\] is 



k = k{t) = -7 



(72) 



which is again a bivalued function. Now, however, the function (p7| ) has no extrema, and it 
is easy to see that we can take a single inverse everywhere as. 



kit) 



-7 



dk 



dt 



so that (|66|), for n = —1, can also be integrated in terms of a zeroth order Bessel function 
i3[ as. 
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Finally putting together the results (^) and (|73|) , we can write 

T_i = 2Ko(ze/3), (74) 

with the parameter (3 given by, 

f^' = A7-^ = A (Ag - A?) . (75) 



Recall however that we have only considered the case A > 0. The case A < involves more 
delicate complex contour integrations, but the final result is the same. 

The integrals X„ with n < —2 are finite in the limit e ^ and they can also be recursively 
calculated from X_i and the following recursion relation for n < —1: 

din n .eAo , . 

^ = -X„-2-.^X„-i. (76) 

This allows us to calculate all the X„ with n < —1 only if we know at least another integral 
besides X_i. Fortunately it is not very difficult to evaluate X_2 from the recursion relation 

^ i din-i 



and by a simple integration. 



e dAr 



A, 



X_2(Ao) = -te / rfA'oX_i(A'o) +X_2(0), 
Jo 

where the integration constant X_2(0), which corresponds to the value of X_2 at the special 
point Ao = 0, can be calculated directly from ( |66D |^ as 

^-2(0)|,>o = ^e-^'^'^^S ^-m\x<o = -JZ^ sin [e(-A)^/2 Ai 



C Geodesic coefficients and xl^ 

We start with either a timelike or spacelike geodesic connecting the points x', x and x, which 
can be written in terms of the proper geodesic distance r as, 

x^' = x'^ (f + e) , x^' = xf" (f ) , x'f" = (f - e) , 

and define at the midpoint a; the parameters t* = d'^t* /dr'^l- and = d"'x^ /dr'^l-, which 
determine the geodesic. From our approximate metric (|2^) in the interaction region, we 
have 

1 f f f ~^ f2 ~^ f2 ~^ f2^ 

J1J2J3 h Jl 1 2 
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and since the metric (|2J) does not depend on coordinates x, y or z, then = Qxx 
Py = gyyUi and Pz = Qzz Zii are constants. Therefore we can rewrite the above expression as 



7-*^ = ^ I I P^ _|_ Py — T?(+*\ (77] 

' /1/2/3 (/1/2P (/2/3)2 (/l/s)^" ^ ^ 

where for convenience we have introduced the function F{t*). Then in terms of and 
derivatives of the function F, the t* coefficients can be written as, 

t; = If, t; = hiF", ti = \ [{ff")' + ff'"] , 
t; = ^ti [{FF"y + FF'"]' , 

where we denote /' = df/dt*. 

Analogously, using that pi = ga x\ are constants, the coefficients can be easily written 
in terms of x\, tl, F and the metric coefficients as. 



x^-^tV x^ ^ jFcf'y + Fg-' 
9 ^ 



{Fg'^'y + Fg''" 



9' 



-t^Xi, 



X. 



1 x\ 
4^ 



All these expressions, evaluated at the midpoint t = t, produce the geodesic coefficients we 
need. 
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Figure captions 



Fig. 1 The colliding wave space-time consists of two approaching waves, regions II and III, 
in a flat background, region IV, and an interaction region, region I. The two waves move 
in the direction of two null coordinates u and v. The four space-time regions are separated 
by the two null wave fronts u = and v = 0. The boundary between regions I and II is 
{0 < -u < 7r/2, V = 0}, the boundary between regions I and III is {m = 0, < f < vr/2}, 
and the boundary of regions II and III with region IV is < 0, v = 0} U {u = 0, v < 0}. 
Region I meets region IV only at the surface u — v — 0. The KiUing-Cauchy horizon in the 
region I corresponds to the hypersurface u + v — Tr/2 and plane wave regions II and III meet 
such a Killing- Cauchy horizon only at V — {u — 7r/2, v = 0} and V' = {u — 0, v — 7r/2} 
respectively. The hypersurfaces u = tt/2 in region II and v = 7r/2 in region III are a 
type of topological singularities commonly referred as folding singularities and they must be 
identifled with V and V' respectively. 

Fig. 2 The subset of Cauchy data which affects the evolution of the quantum field along the 
center u = v oi the plane wave collision lies on the segments {0 < -u < 7r/4, v = 0} U {u = 
0, < V < 7r/4}. Region S is the causal future of this Cauchy data (or equivalently, the 
causal past of the collision center) . 

Fig. 3 We change the mode propagation problem for the plane wave collision, in the causal 
past of the collision center, region 5, by a much simpler Schrodinger-type problem which 
consists in eliminating the dependence of the field equation on coordinate z by taking ^ = 0, 
and substituting the Cauchy data which come from the single plane wave regions on segments 
{0 < < 7r/4, i> = 0} U {rt = 0, < < 7r/4} by much simpler Minkowski Cauchy data. 
This procedure is essentially equivalent to modifying the space-time geometry in the causal 
past of the collision center by eliminating the dependence on coordinate z in the line element, 
setting ^ = 0, and smoothly matching this line element, through plane wave regions II and 
III, with the flat spacetime below the segment {t = 0, — 7r/4 <v< tt/A}. 
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